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ABSTRACT 

This talk presents an introduction to tlie use of dispersion relations to con- 
strain the shapes of hadronic form factors consistent with QCD. The applications 
described include methods for studying \Vcb\ and \Vub\-, the strange quark mass, 
and the pion charge radius. 



1. Introduction and History 

Between the mid 1950s and late 1960s, a great deal of theoretical activity focused 
on attempting to solve (or at least severely constrain) problems of strong interaction 
physics using dispersion theory. An extensive and elegant body of work was developed 
to study the analyticity properties of form factors and scattering amplitudes. Eventu- 
ally, however, when theorists believed they had reached the the limits of what could 
be gleaned from dispersive techniques, their attentions were drawn elsewhere: to the 
quark-parton model, to current algebra, and eventually to gauge theories, especially 
QCD. 

The appeal of dispersion theory lies in its ability to incorporate in a completely 
rigorous and model-independent fashion those features shared by all well-defined field 
theories, namely, causality, unitarity, and crossing symmetry. Moreover, it works equally 
well in perturbative and nonperturbative regimes of the underlying dynamical theory. 
However, no specific Lagrangian is demanded by this scheme, and this lack of speci- 
ficity acts as a double-edged sword: Without dynamical input, one can only deduce 
those consequences common to all possible dynamics. On the other hand, we now pos- 
sess QCD, which is the fundamental, albeit unsolved, theory of strong interactions. 
A combination of the two, in which QCD inputs are inserted directly into dispersion 
relations, should yield a rich harvest of rigorous and model-independent bounds on 
hadronic quantities. In this talk we explore the implementation of this idea to the 
specific cases of weak and electromagnetic hadronic form factors. 

Dispersion theory has been with us in particle physics for quite some time. The 
origin of the name traces directly back to the famous Kramers-Kronig relation^ in 
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electromagnetism, used to describe the dispersion of light in an arbitrary medium. The 
standard formula reads 



where /, uj, a^ot, and P denote, respectively, the forward scattering amplitude, fre- 
quency, total cross section, and principal value prescription to remove the denominator 
singularity. It follows that the dispersive (Re/) and absorptive (Im/ ~ (Jtot) ampli- 
tudes are intimately connected. This relation is derived by using causality and unitarity, 
which lead to restrictions on the analyticity properties of / in the space of complex 
UJ. From there, the elegant theorems of complex analysis, especially Cauchy's theorem, 
provide the identities known as dispersion relations. 

Since quantum mechanics and quantum field theory are expressed over the field of 
complex numbers, it is natural to expect that some variant of the dispersive approach 
should also exist in particle theory. Indeed, as early as 1951, Gell-Mann, Goldberger, 
and Thirringp described how causality and unitarity lead to a dispersion relation for the 
vacuum polarization two-point function of QED. A flurry of other dispersion relations 
followed in the literature, each presented with more or less rigor, depending upon 
assumptions about the analytic structure of the quantity under scrutiny; however, the 
particular dispersion relation used below is nothing more than the QCD version of the 
one first studied in 1951. 

2. Formalism 

We begin by defining the vacuum polarization tensor as the two-point current 
correlator in momentum space: 



Here J is some chosen current; since we are working with QCD, we choose it to be a 
quark bilinear. Moreover, we choose J to be a weak or electromagnetic, rather than 
gluonic, current, so that the individual (perturbative) current insertions are easier to 
identify. Suppressing for now the Lorentz indices ^u, we would like to use Cauchy's 
theorem to write an expression 



which relates 11 at two different momentum arguments, and t. However, in order to 
do this, we must identify a closed contour C, inside of which H is analytic in t. In the 
present case, causality implies that H(t) is analytic in t except on parts of the positive 
t axis, where can create on-shell hadrons, which generates a discontinuity only in 
the imaginary (absorptive) part of H. We choose C to consist of the lower and upper 
sides of this branch cut, together with the circle with \t\ — > oo; the latter contribution 
vanishes as long as H — * for large Then 
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Using the Schwarz reflection principle (11(2;*) = Il*{z) if 11 is real on some segment of 
the real axis, which is true for t < since there are no on-shell thresholds and hence is 
no imaginary part there), the two terms in (^ are equal: 

If n(g^) diverges, or the contribution from the circle \t\ —>■ oo does not vanish, such 
terms may be removed through the process called "subtraction": Since the offending 
terms appear as coefficients of a polynomial in g^, taking a sufficient number n of 
derivatives yields a finite result. 

Then the expression for the dispersion relation reads 

^^^"n! {dq'T ~7tJo (t - g2)n+i ' 

Restoring the Lorentz indices and inserting a complete set of states between J and 
(unitarity) yields 

Imn^'^(t + ^e) = ^E/ d<l>{T){2nr6' (t-E^) (0| Jlr)(r| jt'^lO), (8) 

where only on-shell states F with phase space $ are included in the sum (a consequence 
of reducing the step functions in the time ordering). The matrix elements (0| J'^IF) are 
nothing more than decay constants and form factors — pure hadronic quantities — while 
can be chosen so that n(g^) can be evaluated directly in the fundamental theory 
of QCD. In particular, one chooses to be far from the hadronic (strong coupling) 
region, and then n(g^) may be conujuted using an operator product expansion. A very 
useful observation due to MeimanEl in 1963 is that the jj, = u components of Eq. (^ 
are positive definite, meaning that each hadronic contribution serves only to saturate 
further the partonic (perturbative QCD) side of the dispersion relation. In this way 
one obtains a rigorous inequality between partonic and hadronic physics^ 

One path from Eqs. (0)-(^ leads to the famous QCD sum rules,u which study 
the saturation of the equality between the partonic and hadronic sides. We focus also 
on what this equality tells us about the behavior of matrix elements (0| J'^jF). The first 
work to use the Meiman inequality with QCD inputs was by Bourrely, Machet, and de 
Rafael in 1981. 

As an explicit example, consider the pion electromagnetic form factor: 

{^^{p')\ju^^{p)) = f{<f){p+pr. (9) 

where q = p — p'. Then Eq. (§) becomes 

Imff^(t) >-Lit-Amlf't-'/y{t)\'9it-Aml), (10) 



while the partonic side, finite after two subtractions, is computed to be 



o 



TT 



TT 



+ n.p. 



(11) 



where n.p. stands for nonperturbative corrections such as vacuum condensates. The 
combined inequahty reads 
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In general, one obtains an inequality of the form 

TT Jt+ (t - g2)"+i - ^ 



(12) 
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where t+ is the lowest threshold and Wp is a positive weighting factor arising from 
phase space^and the quantum numbers of the form factor F{t). As discussed by Okubo 
and FushihH in 1971, it is very convenient to map the complex t plane with a cut for 
t+ <t < +00 to the unit disc using a complex kinematic variable z: 



z{t;ts 



(14) 



which maps the^ upper (lower) side of the cut to the lower (upper) half of the unit circle. 
The parameterQ tg < t+ is chosen later for convenience. One then defines a weighting 
function 



<f>F{t;ts 
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\ U(^){q^){t - q^)^+^dz{t;ts)/dt\' 



(15) 



which is analytic inside the unit circle, and in terms of which the dispersive bound 
reads 

:^d-\M^)F{z)f<l. (16) 
Zm J z 

If any subthreshold poles remain inside the unit circle at points z = Zp, they may be 
removed by means of so-called Blaschke factors. 
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1, \Pf(z) 



Pf{z) has the feature that for \z\ 
unchanged, 

^^(j)'^\Mz)PF{z)F{z)\' <l 



1, so that the dispersive bound is 

(18) 



and (j)F{z)PF{z)F{z) is analytic on the whole unit disc. Crossing symmetry relates the 
form factor in all kinematic regimes by analytic continuation. We have thus isolated 
the analytic structure of the form factor, 



where the coefficients a„ are unknown; however, inserting Eq. (|TUp back into (|I^) gives 

oo 

E < 1- (20) 

n=0 

Equations ([19D and (^OD first appearedi in 1995, and re-express in a very compact and 
explicit notation all of the analyticity, unitarity, and explicit QCD information implicit 
in Eqs. and (^. Since (pF and Pf are known functions, the form factor is known 
except for a set of parameters a^, each of which must be less than unity in magnitude. 
A randomly chosen shape for a form factor would almost inevitably have some |a„| > 1, 
and thus would be disallowed by the dispersive bound (pOl) . 

One more point that makes the model-independent parameterization ( ]TP| ) useful 
is that for spacelike and semileptonic processes, the allowed kinematic range for z 
tends to have \z\ -C 1. Indeed, the parameter ts is chosen to enhance this effect. For 
example, for B Dii>, \z\ < 0.03. This means that the convergence of Eq. (|19D is 
geometrically fast, and only the first few relevant to the shape of the form 

factor. The theoretical uncertainty incurred by ignoring the other, infinite set of a„'s 
is called "truncation error, and falls off geometrically fast with the number of a^'s 
used to parameterize the form factor. 

3. A Gallery of Results 

1) \Vub\ o.'iT'd \Vcb\- The^need for a parameterization describing all solutions of 
Eqs. (0)-(|l) was recognized^ in studies of the B iriu form factor, useful for the 
extraction of \Vub\- There it was seen that the inclusion of each (at that time, hypo- 
thetical) form factor data point served to decrease the region allowed by the dispersion 
relation geometrically fast (Fig. |^). Similar comments apply to using points from a 
lattice simulation.El 

The mo del- independent form factor parameterization Eq. ( |l9l) was first useda 
to extrapolate measured B D^*Hu form factor data to a point where phase space 
vanishes. In order to extract | V"cf,| from the form factor^ in the differential width 

(5 ^ D*£u) oc \V,bf F{q') ' J{MB + Mo>y-q\ (21) 



2 ' 



dq 

one must separate \Vcb\ from |F(g^)|. The normalization of F, namely, F{q'^ =AMb — 
Md*Y) = 1 (up to small corrections), is determined by the heavy quark limit.ll3 How- 
ever, Eq. ( ^11) shows that phase space vanishes at exactly this g^; therefore, an extrap- 
olation is needed. Previously, experimental measurements of the form factor used an 

* Strictly speaking, 4 form factors contribute; however, in the liiait of heavy quarks, each one either 
vanishes or is proportional to a single "Isgur-Wise" form factor.ll3 
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Fig. 1. Bounds on the B tt£u form factor / using the dispersive bound and fixing zero, one, 
and two points in (a), (6), and (c), respectively. Dashed hnes indicate pole dominance models. 
(d) shows how certain choices of B* pole parameters can violate the dispersive bounds. 
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CLEO0 
CLEO0 
ALEPH§ 
ALEPH0 
DELPHlllZl 


B D*£v 
B Diu 
B D*(iu 
B Diu 
B D*£iy 


36.9+^:^ 

44.8 ±6.1 

31.9 ±2.4 
29.2 ±7.3 
38.0 ± 1.3 



Table 1. Determinations of \Vcb\ using Eqs. (19) and (po|). Footnotes reference the source of 
the fit. Uncertainties are statistical plus theoretical. 



ad hoc linear or quadratic extrapolation, which implies a theoretical uncertainty of un- 
known ^ze (see Fig. 2a). Using ( [T9| ) and ( pO]) removes this uncertainty, and subsequent 
workll3y refined the analysis to the point that it is used by both theorists and the 
experimental groups themselves (Table p. 

2) Strange quark mass. Ki^ decays possess two form factors, one of which appears 
with the coefficient mf in the rate and is called the scalar form factor. The correspond- 
ing n*^")(g^), evaluated deep in the Euclidean region, is proportional to (m^ — m„)^, i.e., 
is sensitive to m^. One can invert the program of 1), so that a large amount of form 
factor data, thus delineating its shape, is usedti to constrain the function 11'^"^ (g^) and 
hence m<j. Indeed, one finds that F{t) oc m^On, so that (pO]) implies a rigorous lower 
bound on m^. 

Currently, not enough data exists in the world sample for such a determination, 
although DA$NE expects to increase the available pool many times over. Until such 
data exists, one may apply the results of a model, or better, a chiral perturbation theory 





Fig. 2. (a). Fit to CLEO data for the B — > 0*11) form factor using the model-independent 
parameterization (|^) with (sohd) and without (dashed) the dispersive constraint Eq. (20). 
The eye prefers the dashed curve, but it is forbidden by QCD. (6). Fit to timehke pion 
electromagnetic form factor data. Contrast the smoothness of the fit above t = 0.4 GeV^ to 
the oscillations below. 



(xPT) calculation to compute the required shape parameters a„. Believing xPT to a 
level of 5%, 1%, or 1/20% leads to lower bounds on mf^{l GeV) of 40, 90, and 140 
MeV, respectively. 

3) Pion form factor. The parameterization (|I9|) exhibits geometric convergence 
for t < {\z\ < 1). On the other hand, one often possesses data directly on the cut 
t > {\z\ = 1), the timelike region. Does ( |T9|) have anything to say about this region? 
Although one must be much more careful about convergence, the answeili^ appears 
to be yes. Theorems of complex convergence plus knowledge of asymptotic {t — oo) 
properties of form factors allow one to use (^) even in the timelike region. 

For example, for the pion electromagnetic form factor, one obtains the fit of 
Fig. 2b. The presence of the p peak is not put in by hand, but simply emerges from 
fitting to (P^). Note, however, the wild oscillations for t < 0.4 GeV^; one can show that 
these occur due to large gaps in the data for 6 > tt /2 on the unit circle \z\ = 1. These 
large oscillations persist when one analytically continues into the spacelike region, where 
one obtains 

|F(t = 0) I = 2.56 ± 2.00, (r^) = 2.66 ± 3.44 fm^ (22) 

which are rather loose bounds, considering that, e.g., \F{t = 0)| = 1 by charge conser- 
vation. This points to the well-known problem of the instability of analytic continuation 
of discrete timelike data to the spacelike region; now, however one can quantify exactly 
how unstable this continuation is. 

One can also proceed directly with spacelike pion form factor data alone, where 



the geometric convergence of (19) is restored. Then one finds,0 using this model-in- 
dependent parameterization, (r^) = 0.480 ± 0.020 fm^, a few a larger than the usual 
numbers quoted in the literature (~ 0.42 fm^), which use ad hoc parameterizations. 

4. Conclusions 

Dispersive techniques provide an elegant, rigorous bridge between hadronic and el- 
ementary quantities. For semileptonic or electromagnetic decays, they provide a model- 
independent, rapidly convergent parameterization of form factors. We have seen that 
a number of different problems have already been studied using this method. Nucleon 
form factors, the K charge radius, and improvement of timelike form factors are obvious 
future directions. The reader can doubtless imagine many others. 
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